D-dimensions Dirac fermions BEC-BCS crossover thermodynamics 



O 
O 

(N 

> 
O 

m 
> 

o 

oo 
o 

O 



o 

=5 



X 



Ji-sheng Chen* 
Institute of Particle Physics & Physics Department, 
Hua-Zhong Normal University, Wuhan 430079, People's Republic of China 

An effective Proca Lagrangian action is used to address the vector condensation Lorentz violation 
effects on the equation of state of the strongly interacting fermions system. The interior quantum 
fluctuation effects are incorporated as an external field approximation indirectly through a Active 
generalized Thomson Problem counterterm background. The general analytical formulas for the d- 
dimensions thermodynamics are given near the unitary limit region. In the non-relativistic limit for 
d = 3, the universal dimensionless coefficient £ = 4/9 and energy gap A/e/ = 5/18 are reasonably 
consistent with the existed theoretical and experimental results. In the unitary limit for d = 2 
and T = 0, the universal coefficient can even approach the extreme occasion £ = corresponding 
to the infinite effective fermion mass m* = oo which can be mapped to the strongly coupled two- 
dimensions electrons and is quite similar to the three-dimensions Bose- Einstein Condensation of ideal 
boson gas. Instead, for d = 1, the universal coefficient £ is negative, implying the non-existence 
of phase transition from superfluidity to normal state. The solutions manifest the quantum Ising 
universal class characteristic of the strongly coupled unitary fermions gas. 
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I. INTRODUCTION 

Since DeMarco and Jin achieved the Fermi 
degeneracy [1], the property of the strongly coupling 
fermions system under the ultra-cold extreme condi- 
tions serves as a bewitching topic in the fundamental 
Fermi-Dirac statistical physics. 

The universal thermodynamics properties learned from 
the ultra-cold atomic experiments can be related to 
many other realistic many-body physics problems[2, 3]. 
For example, in nuclear physics, the magnitude of the 
neutron-neutron S'-wave scattering length a^vw = —18.6 
fm is much larger than the interaction range given by 
R w l/m-a ~ 1.4 fm, which will indicate the low en- 
ergy thermodynamics universal property in the dilute 
regime [4]. Furthermore, the Feshbach resonance physics 
can be also related with the strongly interacting SU(N C ) 
non-Abelian plasma low viscosity physics realized in the 
ultra-relativistic heavy-ion collisions or existed in the 
core of neutron star [5] and even the "high-temperature" 
T c physics [6]. 

Across the Feshbach resonance regime, the interac- 
tion changes from weakly to strongly attractive according 
to the magnitude of the applied external magnetic field 
B. At this Bose-Einstein Condensation (BEC)/Bardeen- 
Cooper-Schrieffer (BCS) crossover point, the scatter- 
ing length diverges |a| = oo due to the existence of 
a zero-energy bound state. This characteristic leads 
to the universal thermodynamics, i.e., the system de- 
tails do not contribute to the thermodynamics proper- 
ties. The intriguing universal unitary fermions thermo- 
dynamics properties attract much attention recently in 
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the literature [4-11]. 

In the dilute unitary limit, the important physical 
quantities energy density E/N = £ {E/N) ideal and en- 
ergy gap A or critical phase transition temperature T c , 
are the two key parameters to be determined experi- 
mentally and calculated theoretically. The dimensionless 
proportionality constant £ is called the universal factor. 
The A is also proportional to the Fermi kinetic energy. 
Many growing theoretical and experimental efforts about 
the universal factor £ have been made in recent years 
Refs.[4, 7-11] [12]. It is worthy noting that the differ- 
ences for the critical temperature or the energy gap can 
be as large as several times among the theoretical or ex- 
perimental results existed in the literature. Up to now, 
there is not a reliable analytical method to deal with the 
unitary limit problem although there are many updating 
theoretical analytical attempts. Beyond the naive per- 
turbative expansion techniques or the lowest order mean 
field theory, a novel method is to be explored for the 
unitary limit fermions thermodynamics [13]. The essen- 
tial task is how to incorporate the in-medium nonlinear 
quantum fluctuation/correlation effects into the thermo- 
dynamics in a reasonably controllable way. 

Recently, inspired by the homology between this uni- 
tary limit topic and the nonperturbative long range in- 
frared singularity manifested by the hot and dense gauge 
field theory, we make an analytical detour to attack this 
intriguing problem[14]. Our initial observation is that 
the scattering cross section between the two-body par- 
ticles at the zero energy bound unitarity is limited as 
S'-wave a = 47r/fc 2 (with k being the relative wavevec- 
tor of the colliding particles) , while the gauge propaga- 
tor is A M „ ~ g^v/k 2 . The interesting analytical results 
£ = 4/9 and A/ef = 5/18 comparable with other the- 
oretical ones are obtained by taking the still unsolved 
classical Thomson Problem as a potential nonperturba- 
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tive quantum many-body arm. In the non-relativistic 
limit with T = 0, the analytical expression for the en- 
ergy density is the same as that obtained by Steele with 
the power counting in the effective theory framework [1 5]. 

The relativistic continuum quantum field theory for- 
malism provides a natural framework to address the in- 
frared long range singular correlations associated with 
the density of state in many-body physics according to 
quantum statistical physics. In terms of the relativistic 
formalism we can certainly achieve more than those we 
need in a non-relativistic theory, but we can recover all 
the non-relativistic limit physical quantities by expand- 
ing in powers of kf /m at any time especially for the T = 
universal unitary gas thermodynamics. In this work, we 
will continue the relativistic continuum field theory at- 
tempts for the strongly interacting Fermi-Dirac statisti- 
cal physics in terms of the in-medium Lorentz violation. 
After giving the general analytical expressions for the en- 
ergy density functional and pressure as well as entropy 
density according to the density functional theory near 
the unitary limit, the final analytical results compara- 
ble with other theoretical approaches will be obtained by 
applying the corresponding Taylor series expansion. 

Initialed with the recent discussions about the low 
rf-dimensions unitary fermions gas[16], we extend the 
Thomson Problem counterterm analytical method to ar- 
bitrary c?-dimcnsions with the relativistic continuum field 
theory formalism. The d = 2-dimensions solution reflects 
the essential interesting physics. The universal dimen- 
sionless coefficient £ can approach to the specific value 
£ = 0, corresponding to the two-dimensions strongly cou- 
pled electrons with the infinity effective mass m* = oo. 
This is also very similar to the well-known BEC phe- 
nomena for the three-dimensions ideal boson gas with 
the simultaneous vanishing of both the energy density, 
pressure and entropy density at T — 0. Furthermore, the 
energy gap can be zero. Instead, for the one-dimension, 
the £ can be negative, corresponding to non-existence 
of the phase transition from the superfluidity to normal 
state at T = 0. The properties of these solutions man- 
ifest the strongly coupled quantum Ising universal class 
characteristic of the unitary fermions thermodynamics 
realized in the artificial experimental environments. 

The present work is organized as follows. In section 
II, the adopted low energy effective quantum field theory 
formalism of the relativistic continuum Proca Lagrangian 
and corresponding Thomson Problem counterterm non- 
perturbative approach for the strongly interacting uni- 
tary gas thermodynamics are described. The effective 
potential and the general thermodynamics quantities for 
finite scattering length a as well as the concrete com- 
parisons with existed results are presented in section III. 
Moreover, the (i-dimensions universal thermodynamics at 
unitarity are discussed in section IV and the concluding 
remarks are made in the final section V. 

Throughout this work, the natural units with c = H = 
ks = 1 are adopted. 



II. FORMALISM 

A. Effective interaction action 

A nonperturbative approach is crucial to under- 
stand the strongly correlating dense and hot many- 
body physics. Near the strongly interacting unitary 
limit, any form of effective interaction can be used to 
detect the universal thermodynamics properties in the 
pseudo-potential spirit of statistical physics. We ex- 
tend previous work[I4] to near unitary limit regime mo- 
tivated by the successful relativistic quantum many- 
body formalism [1 7]. To model the strongly interacting 
fermions thermodynamics in terms of the fundamental 
in-medium Lorentz violation effects, let us consider the 
relativistic Dirac fermions interacting with an auxiliary 
Proca-like Lorentz massive vector field [18] 

^effective = 4>{il^ - m)V> 

-\-S£ Auxiliary Thomson Counterterm Background- (i) 

The is the vector field with the field stress 

F„ v = d^A v - d v A„. (2) 

In the effective action £/,a-v> = \ m a A ^ A>1 + A* JM > 
the vector current contributed by the fermions is 

r = (3) 

In Eq.(I), m is the bare fermion mass with the vec- 
tor boson mass squared assembling both a bare m 2 A and 
a fluctuating one m^[19], M 2 A = m 2 A + m 2 ^. Here, a 
bare Lorentz violation term l/2rri A A fJ- A f * is included ob- 
viously, which determines the interaction strength in vac- 
uum with the bare electric charge g. For example, the S- 
wave scattering length a is given through the low energy 
Born approximation in the language of the low energy 
quantum mechanics in three-dimensions 

4 = — ■ (4) 

m A m 

B. Induced interaction determined by Thomson 
counterterm background 

The induced Lorentz violation term caused by the 
Thomson background l/2m 2 ^A fjL A 11 represents the many- 
body renormalization effect, which has been introduced 
with a hidden local symmetry manner(yl&eZian Higgs 
model) by mapping the unitary limit topic to the infrared 
singularity correlation in dense and hot gauge theory [14]. 
Taking the coupling constant g to be an "electric" charge 
of the fermions, makes it possible for us to introduce 
the generalized renormalization condition, i.e., the Ac- 
tive Thomson stability condition. With the physical con- 
straint Thomson stability condition, the parameter m 2 ^ 
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is found to be the negative of the gauge invariant De- 
bye mass squared in the unitary limit \a\ = oo[14, 20]. 
The gauge invariant Debye mass associated the density 
of state (DOS) reflects the essential quantum fluctuation 
effect on the thermodynamics. 

More generally the analytical expression for the De- 
bye mass squared is the gauge invariant static infrared 
limit of the 00 ingredient for the vector boson polar- 
ization tensor LT^ (po,p), which can be also calculated 
with the full fermion propagator through the Dyson- 
Schwinger equations or relativistic random phase approx- 
imation (RPA)[21] 
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In Eq.(5), the 0-component of the d + 1-dimensions mo- 
mentum k — (fco,k) in the fermion loop is related to 
temperature T and effective chemical potential fi* via 
ko = (2n + l)nTi + ri* determined by the pole of the full 
fermion propagator with the relativistic Hartree approx- 
imation (RHA) formalism. 

Throughout this paper, the shorthand notation J k = 
J d d k/(2ir) d will be adopted for the d-dimensions ab- 
stract Dirac phenomenology [22, 23]. With finite temper- 
ature field theory, the general analytical expression for 
rf-dimensions unitary fermions gas is derived to be 
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(d-2)25* + k 2 
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In the above expressions, 
1 
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(6) 



(7) 



are the distribution functions for (anti-)particles with 
Ek = Vk 2 + m 2 and (3 = 1/T the inverse temperature. 

This kind of calculations can be indicated by the com- 
prehensive diagrammatic representation as indicated in 
Fig.l. 



C. Diagrammatic expression for the coupled 
generalized Dyson-Schwinger equations 

It should be stressed the in-medium boson propaga- 
tor with the negative of the gauge invariant longitudinal 
L (0, |p| -» 0) = -n^°(0, |p| 0) of the po- 
(0, |p| — > 0) presented in Fig. l.a. 
This minus sign implies that the induced interior corre- 
lation contribution is realized with the external field ap- 
proximation formalism indirectly caused by the opposite 
charged Thomson counterterm background. 

The minus sign in front of LT^(0, |p| — ► 0) is cru- 
cial, which ensures the theoretical thermodynamics self- 
consistency and the unitarity of the final physical re- 
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FIG. 1: The generalized coupled Dyson-Schwinger equations 
incorporating the Landau pole contributions: (a), The full 
fermion propagator determined with the effective interaction 
through the in-medium anti- screened vector boson "propaga- 
tor" in instantaneous approximation oc l/(m\ — n^(0, |p| — > 
0)); (b), The full vector boson propagator calculated with the 
in-medium full fermion propagator through RPA. 



suits in the unitary limit [20]. The comprehensive Lan- 
dau dynamical screening effects are automatically incor- 
porated while the concrete physical contents are quite 
different from the conventional weak coupling resumma- 
tion techniques. With these two Feynman diagrams, 
the relevant algebra equations formalism instead of the 
coupled integral equations one makes it possible for us 
to deal with the low energy density and temperature 
fluctuation/correlation effects on the strongly interact- 
ing fermions thermodynamics in an accurate way. In 
other words, the in-medium density correlation effects 
are included through the effective chemical potential in- 
stead of directly through the effective fermion mass. The 
spirit hidden in this kind of approach is quite similar to 
the Kohn-Sham density functional theory[17, 24]. Es- 
pecially, the integral momentum divergence is naturally 
gauged by the relativistic kinematic factors due to taking 
the density functional theory formalism to approach the 
universal unitary gas thermodynamics. 



That the correlation effects arc embodied through the 
Active opposite electric charged Thomson counterterm 
background can also avoid the theoretical double count- 
ing trouble, which leads to the key minus sign difference 
as found in previous works [20]. This minus sign differ- 
ence is inspired by the anomalous thermodynamics dis- 
cussions of neutron-star nuclear matter with either short 
or long-range force interactions[25, 26]. 
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D. 



Parameters at T - 



At T 
have 



0, from the general expression Eq.(6), one can 



2d 



kjr'E, 



(8) 



which can be reduced to 



k' 



for 



m B _ 2rf "-f 

"s 3- (2v^) d r(|+i) 

the non-relativistic limit. The fermion density in d 
dimensions is 



n = 



h d 

f 



(20F) d r(| + i)' 



(9) 



In above expressions, the T(| + l) is the Gamma function. 



For example, with d = 3, the parameter m 2 ~ is 

JD 



m% 



kfE f 



(10) 



with kf being the Fermi momentum and Ef = k 2 + m? 
the relativistic Fermi energy. In the non-relativistic limit, 
the corresponding result is ^f- = — -^r-. 



For d = 2, one can have 



= — -f- with the ex- 



pression Eq.(6). The corresponding non-relativistic re- 



sult is -gr 



For d 



1, the result is — S- 

' g 2 



2Ej_ 



while the non-relativistic limit is — £■ 



In the 



unitary limit, they are negative of the gauge invariant 
Debye/Thomas-Fermi mass squared[23]. 

Before making the following calculations, let us dis- 
cuss the physical meanings of the coupling constants 
represented by the vector boson mass vha and electric 
charge g. They are the parameters introduced to char- 
acterize the bare vacuum interaction strength (scattering 
length a in three-dimensions) between particles with rel- 
ativistic continuum Dirac formalism. Only in the three- 
dimensions, the coupling constant g can be of dimension- 
less, which can be clearly seen from Eq.(4). Unusually, it 

will have the mass dimensional g', where g' is dimen- 
sionless and \i is an arbitrary mass scale of the problem. 



E. Effective interaction strength 

Physically, the interaction between the fermions is 
rcnormalized by the Lorentz violating many-body envi- 
ronment and the induced interaction strength between 
the particles is characterized by the three parameters vn 2 A , 
m 2 - and q 

B a 



The in-medium effective scattering length 
formalism can give relevant concise expressions to de- 
scribe the strongly correlating effect in a nonperturbative 
manner 

m g 2 1 



2ttM2 
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It is worthy noting that the m 2 ^ and g appear in the 
analytical expressions as a fraction ratio. Actually, the 
magnitude of this fraction ratio m 2 ^ / g 2 characterizes the 
remarkable low energy long range infrared quantum fluc- 
tuating effects associated with the DOS on the univer- 
sal thermodynamics according to the quantum statisti- 
cal physics. Therefore, there are not any remained ad- 
justable or expansion parameters, i.e., the physical vac- 
uum interaction strength parameter/ S- wave scattering 
length a controls completely the final thermodynamics 
quantities as given below. 



III. EFFECTIVE POTENTIAL THROUGH 
LORENTZ VECTOR CONDENSATION WITH 
EXTERNAL FIELD APPROXIMATION 

In methodology, the quantum many-body correla- 
tion/DOS effects are incorporated as an external field 
approximation to ensure the theoretical thermodynam- 
ics self-consistency, from which the effective potential 
can be obtained with the standard path integral by tak- 
ing the relativistic Hartrce instantaneous approximation. 
The external field approximation formalism implies that 
the vector current conservation/gauge invariance can be 
guaranteed by the Lorentz transversality condition real- 
ized by RHA 



M" = 0, 

from which the effective potential reads[17, 27] 
tl/V = -TlnZ 



(12) 



-M 2 A 2 

-M a Aq 



l + e 



-P(E kti -n*) 



+ ln(l + e-^.*+^) 



M\A 2 -2T / [ln(l + e- /3(Bfc - M * ) ) 



+ ln(l + e -^(^+M*))l ; 



(13) 



where "2" represents the hyperfine-spin degenerate factor 
of the fermions system. Furthermore, the self-consistency 
condition or the tadpole diagram with the boson self- 
energy for the full fermion propagator leads to 



A 



g 2-rra, 
M\ 



gm 



(14) 



from which the effective (local) chemical potential /j,* is 
defined with a gauge invariant manner 



2 r Ka e ff 

/i = /i + \xj = ii —n, 



m 



(15) 



where fi is the global chemical potential. The fermions 
particle number (electric charge number) density is 



n 



= 2/ [f-7] 



(16) 
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A. Thermodynamics near the strongly interacting 
unitary limit regime 

With the external field counterterm instantaneous ap- 
proximation formalism of the interior correlating effects, 
one obtains the general analytical expressions for the 
energy density functional from the effective potential 
Eq.(13) with the thermodynamics relation e = y + /j,n + 



— -^-n 2 + 2 

m 



/ E k [f + f] 
Jk 



and pressure P = — 0,/V 



to 



(17) 



(18) 



with the entropy density 

S_ 

V 



= -2 [ [f In / + (1 - /) In [1 - /] + (/ - /)] 

Jk 

K 2 / a(/+/,+ ^l (/+/v "*")' 



(19) 

In Eq.(17) and Eq.(18), the first terms are directly re- 
lated to the interaction and/or quantum fluctuation con- 
tributions. In the weak coupling limit, the bare vacuum 
interaction strength/scattering length a e // — + a can re- 
cover the lowest order conventional mean field theory 
which neglects the quantum fluctuation contributions. 
The second terms in Eq.(17) and Eq.(18) as well as the 
entropy density Eq.(19) appear as very much the ana- 
lytical formalisms for the free Fermi-Dirac gas. However, 
the correlating effects are also implicitly included through 
the effective chemical potential esp. for T ^ 0. 

The Eq.(17) and Eq.(18) with m\ = -11^(0, |p| -» 
0) = n^°(0, |p| — > 0) are the final results. The remaining 
task is to give the numerical results for given tempera- 
ture T and density n with the fermion mass to. In this 
work, we limit ourselves to the T — universal thermo- 
dynamics attempting to obtain the analytical results. 



B. Thermodynamics universality at T = for 
three-dimensions 



At T — and from the above general analytical expres- 
sions, one can obtain the energy density and pressure for 
finite scattering length a 



4 It, k i+ E f 



e = 



TO 



2k f E;a 



TO — 



2k, E 



-n z + 



{2k 2 +m 2 )k f E f -m i In 

8^2 ' 

(2k 2 - 3m 2 )k f E f + 3to 4 In 



24tt 2 



(20) 



In the non-relativistic limit with T — 0, the analytical ex- 
pression for the energy density can recover Steele's main 
result obtained within the effective theory framework [15]. 
In deed, the non-relativistic negative Debye/Thomas- 
Fermi mass squared can readily approach the non- 
relativistic power counting result of Steele with the gen- 
eralized coupled Dyson-Schwinger equations (through the 
full anti-screened vector boson propagator) as indicated 
by Fig.l. The non-trivial physics occurs at the pole of 
the first term in Eq.(20), i.e., in the repulsive BEC regime 
with to — 2kfEfa/n = — > 2kfa/-K w 1 as pointed out 
in Refs. [13, 15]. 

Due to the universal properties at T = in the dilute 
unitary limit, the energy density can be scaled as [9] 



5 f I ? k f a 3 k 2 a 2 



0( 



(kfa) 



(21) 



where £f = k 2 /(2m) is the Fermi kinetic energy. Apply- 
ing the Taylor series expansion of e Eq. (20) according to 
l/(akf) as well as kf by keeping only up to the lowest 
order of kf/m, the corresponding universal coefficients 
£ = 4/9 = 0.44, C = 5tt/18 = 0.87, v = n 2 /l2 = 0.82 are 
exactly consistent with those £ w 0.44,^ w l,v w 1 of 
Refs. [7, 9] obtained with Monte Carlo calculation. Espe- 
cially, one of them with much attention is £ = 4/9. The 
recent quantum Monte Carlo calculation also gives the 
same result [28]. 

In the unitary limit, the strong fluctuation and corre- 
lation effects not only do manifest themselves in the bulk 
properties but also modify the quasi-particle properties 
in a substantial way from the viewpoint of the univer- 
sal thermodynamics. At unitary \a\ = oo, the energy 
gap A is derived from the total energy density with the 
odd-even staggering (OES) 



A = / s 



2kj + 3to 2 - 3mEf 
' 6^ ; 



(22) 



with f sw being the Fermi-Dirac statistical weight factor. 
In the non-relativistic limit with the factor f sw = 5/3, 
the analytical result is A/e/ = 5/18, which is reasonably 

consistent with theoretical result -4-e/e 2k f a ~ 0.40e/ 
obtained through the BCS model but with the effective 
scattering length. The corresponding phase transition 
critical temperature is still approximated by the BCS re- 
lation T c — e 7 A/-7r with 7 being the Euler-Gamma con- 
stant 



T c w 0.157T/, 



(23) 



as given in Ref.[14]. This is in reasonable agreement with 
the updating results[10, 28-30]. 

The pressure is P = 1/QPfg m the unitary limit, from 
which one can find the sound speed is reduced remarkably 



v = 



V F G = 



V2 



-Vf, 



(24) 
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where vp G = l/3v 2 is the sound speed squared for the 
ideal Fermi gas with the Fermi velocity Vf = kf/m. 
For the ideal fermion gas, the ratio of pressure to en- 
ergy density is well known 2/3. For unitary fermion gas, 
this ratio is found to be changed to 1/4 for the non- 
relativistic occasion due to the strongly correlating ef- 
fects. This is different from that obtained with the uni- 
versal hypothesis based on assuming the scaling property 
in terms of the zero-energy bound state. This difference 
can be attributed to the implicit pairing correlation con- 
tribution to the binding energy in the strongly coupling 
limit [31]. The simultaneous experimental detection for 
the universal coefficient £ and the sound speed can judge 
this dilemma. We also note that similar conclusion can 
be seen in the analytical Dyson-Schwinger attempt of 
Ref.[30]. 



IV. D-DIMENSIONS FERMIONS UNIVERSAL 
PROPERTIES AT UNITARITY 

Although the unitary limit issue is proposed for the 
three-dimensions physics representing for such as the in- 
ner neutron star crust physics in the low energy strongly 
interacting nuclear many-body theory context, it is very 
interesting to visit the d-dimensions occasion of the gen- 
eral Fermi-Dirac statistics as a toy model in the unitary 
limit by applying the Thomson Problem counterterm an- 
alytical method. 

Concentrating on the realistic physics realized in the 
experimental environments, the non-rclativistic analyti- 
cal results reduced from the general expressions Eq.(17) 
and Eq.(18) in the unitary limit with the Thomson Prob- 
lem counterterm approach are given below. The ki- 
netic energy density contribution for d-dimensions non- 
relativistic fermion gas is 



gap 



- i<l( 



2d 



e f kj 



(2VSFf (d + 2) r(f + i)' 



while the pressure is 



Pideal 



Sfk 



(2^) d (d+2)ni + iy 

The negative correlation energy is derived to be 

2 e/fc/ 



£int 



(2y^) d dr(l + 1) 



(25) 



(26) 



(27) 



Analogously to the discussion for three-dimensions, 
one can obtain the general universal dimensionless co- 
efficient £ expression according to dimension d 

(d-2)(rf+l) (rf-2)(d+l) 

? - J sw JU I o\ - ^2 • ^1 



d(d + 2) 



With the OES method for the d-dimensions unitary 
fcrmions gas, one can have the extended S-wave energy 



A/e, = f s 



d-2 _ (d-2)(d + 2) 



2d 



2d 2 



(29) 



In above expressions, f sw = (d + 2) jd is the d- 
dimensions Fermi-Dirac statistical weight factor in the 
non-relativistic limit. It is worthy noting that the ratio 
of pressure to energy density is changed from 2/d of the 
ideal fermion gas to 1/(1 + d) of that in the unitary limit. 

From Eq.(28), one can find that £ will be equal to zero 
at d = 2. This corresponds to the infinity m* = oo of the 
effective fermion mass. At this specific two-dimensions, 
the pressure as well as entropy density is also equal to 
zero, which is quite similar to the BEC phenomena for 
the ideal boson gas at three-dimensions. In this specific 
occasion, the sound speed will approach to zero in terms 
of the Landau Fermi liquid theory. Furthermore, the en- 
ergy gap will be zero in the meantime with the OES result 
Eq.(29) of the two-dimensions. This result indicates the 
low two-dimensions unitary limit thermodynamics prop- 
erties are exotic. 

In the strongly coupling limit, the essential charac- 
teristic to* — ► oo is in line with the strongly coupled 
two-dimensions electrons discussed in the literature [32]. 
The zero energy gap is so similar to the vanishing of 
the extended S pairing correlation etc. with the increase 
of repulsive interaction strength in the two-dimensions 
Hubbard Model with quantum Monte Carlo study [33]. 
It is worthy noting that the behavior of £ according to 
the spatial dimensions d as indicated by Fig. 2 is different 
from those found in the recent works[16]. 




FIG. 2: The universal coefficient £ versus dimension d. The 
lower solid line (a) is for the non-relativistic Eq.(28), while the 
upper dotted line (b) is for the ultra-relativistic limit Eq.(30). 
The dashed parts correspond to the instability region. 

In Fig. 2, we also give the ultra-relativistic limit result 
for completeness with the explicit analytical expressions 

(d-l)(2d+l) (d-l)(2d+l) 



£' = f 

S J SV 



2d(d+l) 



2d 2 



and 



A'/4 = K, 



(d-1) (d-l)(d+l) 



2d 



2d 2 



(30) 



(31) 
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It is worthy noting that in this opposite limit the cor- 
responding correction of the Fermi kinetic energy is the 
magnitude of the Fermi momentum kf, i.e., e'j — > kf. 
The prime symbol has been adopted in order to avoid the 
confusion with the non-relativistic results. In Eqs.(30)- 
(31), the statistical weight factor is f' sw = (d+l)/d. 
There exists a rummy duality relation between the non- 
relativistic and the ultra-relativistic results seen from 
above expressions. 

From Fig. 2, one can clearly see that these solutions 
manifest the quantum Ising universal class characteristic 
of the strongly interacting fermions gas realized in the 
experimental environments [34]. For the non-relativistic 
occasion, the universal coefficient is £ < for d < 2 with 
the negative sound speed squared. The negative sound 
speed squared means that there is no phase transition 
from the superfluidity state to normal state due to the 
spinodal instability of the superfluidity state. In the one- 
dimension, the strong fluctuation effects in the system 
will make the regions out of phase with each other and 
reduce the possibility of reaching the thermodynamics 
equilibrium. This is consistent with the general phase 
transition theory, i.e., there is not a phase transition for 
one-dimension. 

With the analytical results such as Eq.(28) and Eq.(29) 
as well as Eq.(30) and Eq.(31), one can also see that the 
proportionality coefficients do depend only on the spatial 
dimensions d and are independent of the system details in 
the unitary limit. This is of the genuine thermodynam- 
ics universality characteristic associated with the critical 
phenomena in terms of the general statistical physics. 

V. SUMMARY 

In conclusion, we have given the general analytical for- 
mulas of the energy density functional and pressure as 
well as entropy density at both finite temperature and 
density in terms of the relativistic quantum many-body 
formalism near the unitary limit regime. The interior 
non-linear fluctuation/correlation effects on the thermo- 
dynamics have been incorporated with an external field 
approximation scheme through a fictive opposite charged 
Thomson background. This detour will allow us to go 
beyond the lowest order mean field theory or the naive 



pcrturbative expansions in an accurate way. 

The d-dimensions universal thermodynamics in the 
unitary limit are discussed. It is found the thermody- 
namics properties of low dimension d — 2 are very in- 
teresting in the unitary limit, i.e., the universal dimcn- 
sionless coefficient £ = with simultaneous vanishing of 
energy density, pressure, energy gap and sound speed as 
well as entropy density. In the non-relativistic limit, the 
effective fermion mass can approach to m* = oo, corre- 
sponding to the strongly coupled electrons in condensed 
matter physics. This characteristic is quite similar to 
the well-known BEC of the ideal boson gas at T = for 
three-dimensions. These d-dimensions solutions mani- 
fest the quantum Ising universal class characteristic of 
the strongly interacting low temperature fermions gas 
physics. 

In addition to the exact consistency of the obtained 
physical results with some existed theoretical analyt- 
ical or simulation ones, this attempt with the un- 
known side/Thomson Problem as a potential quantum 
many-body nonperturbative arm to solve the other un- 
known side facilitates the concrete comparison of the 
non-relativistic many-body methods and the Dirac phe- 
nomenology of the relativistic continuum formalism in 
terms of the fundamental in-medium Lorentz violation. 

The essential features of the strongly coupled unitary 
gases may be captured by the Lorentz violation with vec- 
tor condensation formalism. The key physics may be ob- 
tained by the external field realization scheme of the in- 
terior correlations between the strongly interacting parti- 
cles through a fictive Thomson background but still with 
the standard RHA and RPA techniques (or the general- 
ized coupled Dyson-Schwinger equations). 
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